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ABSTRACT 

In studies of star-forming regions, near-infrared excess (NIRX) sources — objects with intrinsic colors 
redder than normal stars — constitute both signal (young stars) and noise (e.g. background galaxies). 
We hunt down (identify) galaxies using near-infrared observations in the Perseus star-forming region 
by combining structural information, colors, and number density estimates. Galaxies at moderate 
rcdshifts (z = 0.1 - 0.5) have colors similar to young stellar objects (YSOs) at both near- and mid- 
infrared (e.g. Spitzer) wavelengths, which limits our ability to identify YSOs from colors alone. 
Structural information from high-quality near-infrared observations allows us to better separate YSOs 
from galaxies, rejecting 2/5 of the YSO candidates identified from Spitzer observations of our regions 
and potentially extending the YSO luminosity function below K of 15 magnitudes where galaxy 
contamination dominates. Once they are identified we use galaxies as valuable extra signal for making 
extinction maps of molecular clouds. Our new iterative procedure: the Galaxies Near Infrared Color 
Excess method Revisited (GNICER), uses the mean colors of galaxies as a function of magnitude to 
include them in extinction maps in an unbiased way. GNICER increases the number of background 
sources used to probe the structure of a cloud, decreasing the noise and increasing the resolution of 
extinction maps made far from the galactic plane. 

Subject headings: dust, extinction — ISM: structure — galaxies: colors — stars: pre-main sequence 



1. INTRODUCTION 

The near-infrared (NIR) has been a valuable window 
into star- forming regions, providing us with a variety of 
tools to study the process by which dark clouds coalesce 
into stars: luminosity functions in the K band (2.2/mi) 
can provi de an accurate estimate of the initial m ass func- 
tion (e.g. lMuench et al.ll2002HStolte et al.ll2005D : studies 
of embedded clusters reveal important details of the pro- 
cesses by which stars form either in groups or in isola - 
tion (e.g. lLada et al.lll99l1 : iRoman-Zuniga et aTJ I2007T ): 
and extinction mapping in the near-infrared allows pre- 
cise deter mination of the column density structure of the 
cloud fe.g. lAlves et al.ll2001t ICambresv et al.ll2002D . 

Of particular use is the narrow range of intrinsic near- 
infrared colors of main-sequence stars (typically H — K 
= - 0.4 and J — H = - 1.0). This narrow range is 
due both to the near-infrared bands lying on the Wien 
tail of all stellar (hydrogen-fusing)-temperature black- 
bodies and to the relative paucity and uniformity of ab- 
sorption feature 1 . By assuming that all stars have the 
same intrinsic color we can do purely photometric ex- 
tinction mapping. This is the heart of the Near Infrared 
Color Excess method (NICE) (|Lada et al.lll994h and the 
Near Infrared Color Exce ss method Revisited (NICER) 
ijLombardi fc Alvesl[200lh . 

Principally H — , but also CO features which are sensitive to 
surface gravity and produce a small split between dwarf and giant 
colors 



Additionally, we can use the narrow range of stellar 
near-infrared colors to identify young stars by their near- 
infrared excess. The thermal contribution from their disk 
or envelope changes their color from that of a plain pho- 
tosphere and places them in a certain region of a near- 
infrared color-color (J — H ve rsus H — K) diagram . This 
is the CTTS locus defined bv lMever et all (|1997l ). 

Unfortunately, a number of other astronomical objects 
also have intrinsic red colors, similar to CTTS. We re- 
fer to all intrinsically red objects as near-infrared ex- 
cess (NIRX) sources. Of particular concern for stud- 
ies of star forming regions are objects which are not ei- 
ther young-stellar objects (YSOs) or T-Tauri type young 
stars (TTS). These objects — brown dwarfs, a variety of 
evolved stars, galaxies, and AGN, are often studied in the 
near-IR in their own right, but in studies of star forming 
regions where we have limited information about their 
nature (perhaps only J,H, and K photometry) they are 
contaminants which must be understood and identified 
in order to avoid common biases in these studies. For 
example, including intrinsically red objects in NICER 
produces an overestimate of the extinction. 

In this paper we identify as galaxies a large number 
of NIRX sources at moderate redshifts (z = 0.1-0.5) in 
high-quality near-infrared images outside of the Perseus 
molecular cloud complex. After describing our observa- 
tions and reduction in fj2]we present a detailed case that 
our NIRX sources are galaxies in In the heart of the 
paper we show how deep NIR images can help identify 
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genuine YSOs in Spitzer observations (2J and how we 
can make use of these contaminating galaxies as addi- 
tional probes for extinction mapping (SJ5J). 

2. OBSERVATIONS 

2.1. Acquisition and Reduction 

The observations discussed herein are deep near- 
infrared observations of the Perseus Molecular Cloud 
complex, obtained as part of the COMPLETE (Co- 
ordinated Molecular P robe Line Extinct ion and Ther- 
mal Emissi on) Survey (iRidge et al.l I2006D and were in- 
troduced in lFoster fc Goodmanl ( 20061 ). These observa- 
tions were obtained at the Calar Alto Observatory on 
the 3.5 meter, using the OMEGA 2000 wide-field NIR 
camera. We observed 6 fields in the cloud (four covering 
B5 and one on each of two dark clouds in the south-east: 
L1448 and L1451), as well as two control fields outside of 
the boundaries of the complex in regions of low extinc- 
tion as estimated from a 2MASS-based extinction map 
(|Ridge et al.ll2006l ). OMEGA 2000 offers a 0.45"pixel 
scale and low image distortion over the full 15.1'xl5.1' 
field of view. Typical seeing conditions were ~ 0.6-0.7". 
We observed in three bands with short, dithered expo- 
sures. Total integration times per field were 45 minutes 
at H , 37 minutes at K s and 6 minutes at J. The control 
fields were observed with the same exposure times as the 
on-cloud frames, and reached roughly the same depths 
(J, H ~ 20 and K s ~ 19) 

Following standard instrumental flat-field and dark 
corrections, images were combined using the xdimsum 
package under IRAF. This package performed sky sub- 
traction, cosmic-ray removal, and image co-addition. 
With our dither pattern, the end result is a trimmed 
image with roughly constant signal to noise over a 13' by 
13' region. 

Photometry was performed on these trimmed im- 
ages in two different modes. Source Extractor 
(jBertin fc Arnoutslll996[ ) was run on the images with a 
5cr detection threshold and an estimated FWHM of 1.15" 
input (see i)3.1.4i Figure[6]) to make a distinction between 
extended and point-like objects. Along with photometry 
in the three observed bands, a classification parameter, 
C, is returned by Source Extractor. C is a number be- 
tween and 1, with 1 being a perfect point source. Our 
basic classification separates objects into point like or 
extended depending on their value of C being larger or 
equal or smaller than 0.5. 

Source Extractor photometry produced a small num- 
ber of very strangely colored (i.e. off the boundaries of 
the color-color diagram shown in Figure[T]either blueware 
or redward) objects which, on inspection, were objects 
at the edge of the frame with an unreliable sky estimate, 
stars in the wings of saturated bright stars, un-cleaned 
cosmic rays or bad pixels, or pieces of the diffraction 
spikes of bright stars misidentified as genuine objects. 
To eliminate such objec ts, we used the Phot Vis package 
(|Gutermuth et al.ll2004D which performs DAOPHOT de- 
tection and aperture photometry with a fixed sky annu- 
lus and aperture size (here 5 pixels for a typical stellar 
FWHM of 2.5 pixels). The Phot Vis package allows easy 
display and manual addition or removal of objects, so a 
hand-cleaned catalog was produced for each image. 

For the control fields, the catalogs from both programs 
were cross-referenced to find objects with detections in all 



three (J,H, and K s ) Phot Vis catalogs and a detection in 
the i?-band Source Extractor catalog. i7-band was cho- 
sen for the classification image as it was the deepest and 
hence highest signal to noise image. For the data fields, 
the requirement of a detection in J-band was dropped, 
as differential extinction within the cloud caused this 
band to drop out for many sources. Magnitudes dis- 
cussed herein are aperture photometry magnitudes from 
Phot Vis, and are not necessarily accurate for extremely 
bright stars H mag at H) or very extended objects 
(> 2.5 " at H). Our control fields contain only a few 
such objects, ~ 10 total or 20 D -1 . Photometric er- 
rors are estimated from the sky noise in PhotVis. The 
control-field catalogs analyzed here were trimmed to in- 
clude only sources with high enough quality colors to ac- 
curately gauge their position on the color-color diagram. 
This definition was o~ C oior < 0.15 mag where 
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Images were calibrated using 2MASS stars within the 
field. This produces a convenient way to establish a pho- 
tometric system without the need for standard star ob- 
servations and airmass corrections. The 2MASS point- 
source catalog (jSkrutskie et al.ll2006h was queried for all 
objects with high quality photometry (flag = AAA), 
cross-referenced with detected objects and used to de- 
termine a coordinate system. 2MASS stars with magni- 
tudes between 11.5 and 13.5 were then linearly fit against 
detected magnitudes to establish a conversion between 
instrumental and 2MASS magnitudes. The bright cut- 
off was chosen to eliminate potentially saturated stars or 
ones where the central pixel counts indicated the possi- 
bility of non-linear detector response. The formal uncer- 
tainty of this fit is not incorporated into the estimated 
errors for these stars, but constitutes an additional sys- 
tematic error term. The H and K s filters on OMEGA 
2000 are the 2MASS filters. The J-band filter is some- 
what broader at long wavelengths (a 50% cutoff at 1.345 
/zm compared to 1.4 ^m for the 2MASS filter), but is 
otherwise similar. Because of the similarly of filters and 
our calibration on 2MASS objects, we consider our ob- 
servations to be on the 2MASS system and do not apply 
additional color-transformations. 

2.2. A clump of NIRX sources 

The essential puzzle we address is illustrated in Fig- 
ure [TJ in both control fields there is a set of objects with 
fundamentally non-stellar colors - colors which cannot 
be produced by any amount of reddening of a star. We 
define Near-Infrared Excess Sources (NIRX) in the fol- 
lowing way: we offset from the tip of the main-sequence 
stellar locus (spectral type M6) by a buffer comparable 
to our photometric uncertainty and draw a reddening 
vector through this point. All objects to the right of this 
line are considered NIRX sources (see shaded region in 
Figure [I]). There is a distinct clump of objects in the 
same location in both of our control fields. This set of 
objects coincides roughly with the region of the diagram 
where young stellar objects (YSO) lie. Given that Fig- 
ure [1] shows control fields far from the cloud that should 
not be littered with YSOs, what is the true nature of 
these objects? 
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Fig. 1. — Near-infrared color-color diagram for (a) Control 04 and (b) Control 01 in Pers eus. The stel l ar loc us is shown, consisting of an 
upper (giant) bra nch and lower (dwarf) b ranch, along with the classical T-Tauri locus from Meyer ct al. (1997). Also displayed is a typical 
reddening vector (Rickc & Lcbofsky 1985). We identify as Near-Infrared Excess (NIRX) sources all objects falling to the left of a reddening 
vector drawn through the tip of the dwarf stellar locus. Photometric errors will scatter some objects across this line in both directions. 



Source Extractor classifications are reliable for identi- 
fying bright galaxies with significant extent on the sky. 
This distinction becomes more difficult when there are 
image distortions toward the edge of a frame (a prob- 
lem for many large detectors), in dense or confused 
fields, when galaxies are small or faint, or when there 
is significant nebulosity in the star-forming region. The 
last is likely to be always true for deep enough images 
due t o the presence of "Cloudshine" ([Foster fc Goodman] 
2006), that is, ambient galactic starlight reflecting off 
dust grains. 

Faint sources constitute the bulk of our NIRX sources. 
Figures [2] & [21 show sources from Control 04 and Control 
01 in three different ways. In Figure^, symbol types are 
assigned based on whether C (classification parameter) 
is greater than or less than 0.5. Our estimated complete- 
ness limits (discussed in § I2.3|) are shown. The bulk of 
sources pile up at the extreme values of C. The faint ob- 
jects in the middle (colored green) can not be assigned 
with the same level of confidence. Figure^ shows that 
there is no significant clustering of extended objects, nor 
are they more prevalent at the edges of the field, indi- 
cating that image distortions remain small towards the 
edge of the field. Figure shows that this division of 
extended versus point-like sources largely corresponds to 
a separation between sources near the stellar tracks and 
our clump of near infrared excess sources, though there 
are exceptions. To account for sources with marginal 
classification we consider the structural information from 
Source Extractor in combination with number density es- 
timates of possible NIRX sources to securely identify the 
nature of our objects. 

We analyze in detail just one ("Control 04") of the 
control fields. This field is significantly deeper in K due 
to a combination of rejected frames and poor weather in 



the other ("Control 01" ). In Control 04 we have sufficient 
quality and depth to make a very secure case for the 
nature of the near infrared excess sources. 

2.3. Completeness Estimate 

To estimate the contribution of various objects with 
different luminosity distributions (e.g. YSOs or galaxies) 
to our survey requires knowledge of the completeness of 
our control field catalogs. This is complicated by the re- 
quirement that we detect a source in all three bands, that 
the error on its color be small (< 0.15 mag), and by our 
source detection algorithm, which is a cross matching of 
the automated Source Extractor and aperture photome- 
try from the Photvis package with some hand selection 

Rather than attempt to simulate this, we begin by ex- 
amining the cumulative histogram for both extended and 
point-like sources in K . From the turnover in these plots 
(see Figure [4j we consider our point-like sources to be 
complete to 18th magnitude and our extended sources to 
be complete to 17.4 magnitudes. Our estimated Ay of 
0.6 corresponds to an Ak ~ 0.07, and is small enough to 
be ignored when comparing number counts in our images 
to if-band magnitude distributions of candidate objects. 
We analyze the subset of sources brighter than these com- 
pleteness limits in [J3] 

3. COLOR DESCRIPTION AND NUMBER ESTIMATE OF 
CANDIDATE NIRX SOURCES 

To test the reliability of our structural information we 
seek a priori estimates for the number of NIRX sources 
we expect to see. Various point-like objects can have 
NIRX colors: young stars (already discussed), evolved 
stars (with a shell of cool ejected material), brown dwarfs 
(particularly young L-type dwarfs in the solar neighbor- 
hood) and quasars. Several kinds of sources can have 
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Fig. 2. — Most near-infrared excess sources in Control 04 are extended objects. Panel (a) shows the classification of objects from Source 
Extractor; also shown are estimated completeness limits for extended (open/red) and point-like (closed/blue) object in K, at 17.4 and 18 
magnitudes respectively. We highlight ambiguous classifications in green. Panel (b) shows that neither type of objects is clustered spatially. 
Pancl(c) shows how this division captures the difference between stellar and NiRX sources in a color-color diagram. A typical galaxy/star 
line of J — K = 1.3 is also shown, but would misclassify a number of sources. 



extended profiles and be NIRX sources: normal galaxies 
at z > (normal local galaxies, composed only of stars 
and dust, should have either stellar or reddened-stellar 
colors), or exotic galaxies such as star- forming or AGN- 
dominated galaxies. Finally, visual binaries may or may 
not show up as extended depending on separation, and 
may or may not have stellar colors. 

In Figure [5] we show a, J — H vs H — K color-color 
diagram with a literature compilation of colors of differ- 
ent galactic and extragalactic objects which are potential 
NIRX sources: the bright est (i.e. local) 2MA SS galax- 
ies from the compilation of lJarrett et al.l (|2003l ): galaxies 
with signi ficant star f ormat ion and dust emission from 
th e list oflHunt et alj ((2002) ; L dwarf stars from the list 
of iKoen et al.l ( 2004h : T-tauri stars with different am- 
m ounts of ext i nction and excess from the observations 
of lEiroa et alj (12001D; colors of Sey fert galaxies from the 
list ofE onso-Herrero et all (11998^); an average sequence 



TABLE 1 

Number Estimates for Objects in Control 04 





Data 




Estimated Contributions 


Point-like 


NIRX 


26 


18 BDs, 2 Quasars, YSOs, AGBs 




normal 


818 


Normal Stars 


Extended 


NIRX 


87 


140 galaxies (not all NIRX) 




normal 


17 


4 visual binaries 



of T-d warf colors constructed bv IZapatero Osorio et aTl 
(|2007l ) from literature data; and the regions typically oc- 
cupied by evolved stars: carbon stars and Long-Period 



Variable stars (with p eriod > 350 days) as defined by 
iBessell fc BretH (|1988l ). 

In the following section we estimate the contribution of 
the different kinds of objects listed above to our observed 
field. A summary can be found in Table [TJ 

3.1. Galactic candidates 
3.1.1. YSOs 

If an observed field is known to have recent or cur- 
rent star formation, NIRX so urces are often t a ken t o 
be YSOs. As explained by lLada fc Adarnsl (|1992f) . 
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Fig. 3. — The same breakdown as shown in Figure[2]for Control 01. The same behavior is seen as in Control 04, but the K image is less 
deep, providing fewer sources and larger errors on the H — K color, and thus decreasing the clustering of sources in color-color space. 
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Fig. 4. — Cumulative histogram and estimated completeness of 
point-like (blue) and extended (red) sources in Control 04. 

the combination of the stellar photosphere and disk 
place YSOs as near infrared excess sources, an exam- 
ple of which is a classic T-Tauri locus (CTTS locus) 



of lMever etail (|1997f) . These sources are often used to 
trace the distribution of star formation within molecular 
cloud s and to identify young embedded clusters (iLi et alJ 
| 1997tlRoman-Zuniga et al.ll2007t iGutermuth et alj|2005t 
Ferreira fc Ladall2007i r In our case, while there certainly 
are a number of young stars in the cloud regions, it seems 
unlikely that there is a large population of hitherto un- 
known young stars with disks present in both our control 
field which were selected to be outside the main regions 
of star formation in Perseus. 

Without a significant degree of clustering, the sources 
would not be a new young cluster, but a population of 
escaped sources from nearby star-forming sites, a situ- 
ation that would raise current estimates of typical star 
formation efficiency or dispersal of young stars greatly. 
Assuming a distance of 250 pc to Perseus, Control 04 
is ~ 9 parsecs away from the nearest dense portion of 
Perseus (B5) and ~ 14 pc away from the nearest site of 
clustered star- formation in IC348. I f a YSO were born 
with a typical velocity of 1 km/s (jKroupa fc Bouvierl 
2003), it would take between 9 and 14 Myr to reach 
the control field areas, which is larger than the typi- 
cal peri od of NIR emissio n of a circumstellar disk (~ 
5 Myr (jHaisch et al.ll2001h ). Some small percentage of 
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Fig. 5. — The locations of candidate NIRX sources as taken from the literature. The thick dark line sho ws the stellar locus, w hile th e 
purple-dashed line s hows the classical T-Tauri locus. In (a) we sh ow ca ndi dates within the galaxy: YSOs ( i|3.1.1> . evolved stars ( i|3.1.2t . 
and brown dwarfs ( £|3. 1.31 1. In (b) we show extragalactic sources ( EI3.2.2I &; i|3.2.1l l 

stars may leave their birth sites with significantly higher 
velocity, and so reach our control field before their disks 
evaporate. However, to constitute a significant number, 
the star formation rate in current clusters would have be 
dramatically higher in order that these velocity outliers 
would still be found in any number 14 pc away. 



3.1.2. Evolved Stars 

A star may also become a NIRX source at the end 
of its lifetime as it sheds its outer atmosphere, produc- 
ing a shell of cold dust around the hot central star. 
This is the same basic mechanism by which a YSO be- 
comes a near infrared excess source. This is a broad 
category, encompassing planetary nebulae, carbon stars, 
AGB sta rs, various variable sta rs, and post-AGB stars. 
Recently, lLombardi et al.l (|2006i) suggested that the ma- 
jority of contaminating near infrared excess sources in 
their study of the Pipe Nebula are bright AGB stars, lo- 
cated at a similar distance near the galactic center (as 
evidenced by a narrow range in magnitudes). With both 
our fields observed far from galactic center the contami- 
nation from these objects would be minimal. Addition- 
all y, our objects a r e muc h fainter than those reported 
by lLombardi et al"] (|2006f) . Other evolved stars are rela- 
tively rare, and are unlikely to contribute significantly. 

3.1.3. Brown dwarfs 

The near-IR c olors of brown dwarf s are p roblematic. 
As illustrated bv lStephens fc LeggettJ (|2004[ ). the colors 
of brown dwarfs depend sensitively on the photometric 
system in which they are measured, since the same ab- 
sorption features (H2O) dictate the flux of the brown 
dwarfs, the design of NIR filters used from the ground, 
and the transparency of the atmosphere on a given ob- 
serving night. 



iKnapp et al.1 (|2004h show the co lors of the L and T 
dwarf s in the MKO system, and IStephens fc Leggettl 
(2004) provides the transformation between 2MASS and 
MKO. Within the uncertainties, it seems that L-dwarfs 
lie in the NIRX region of our color-color diagram, while 
T-dwarfs do not (see Figure [5]). 

We estimate the expected num ber of L-dwarfs as fol- 
lows. IChabrier fc B araffe (2000) estimate the galactic 
di sk brown dwarf d ensity at 0.1 pc~ 3 and Figure 12 
of I Cruz et al.l (|2003( ) provides a luminosity function for 
nearby ultracool objects (including M7 - L8) selected in 
a region of color-color space similar to the location of our 
NIRX sources. Even the brightest L-dwarfs are visible 
in our survey [K = 18 mag) only if they are relatively 
nearby. For each magnitude bin we estimate the fraction 
of the total population in this bin, calculate the distance 
out to which this group is observable for our survey, mul- 
tiply the derived volumes by the number density, and 
sum over all magnitude bins. This provides a very con- 
servative estimate of the contamination, as we include 
many M7 - M9 stars whose colors lie partially overlap- 
ping the color locus of non-brown dwarf main-sequence 
stars. Our upper bound estimate is ~ 18 brown dwarfs, 
a significant contribution. 

3.1.4. Binaries 

A large fraction of our stars are expected to belong to 
binary systems. An unresolved binary star may not lie 
on the unreddened stellar locus, but its color should still 
lie generally close to the locus as the average of two or 
more unreddened colors should stay close to the range 
of the separate components. Chance superpositions of 
two stars (i.e. a visual binary) could result in Source 
Extractor classifying such an object as extended. Such a 
superposition would be an extended object with a star- 
like color. 
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In order to estimate the number of such superposi- 
tions we carried out a scries of Monte Carlo Simulations. 
To include the possibility that faint stars would distort 
the profiles of brighter objects, we lowered our detection 
threshold to 3c in H, producing ten times as many total 
sources as we analyzed in the rest of this paper. We mea- 
sured the point-spread function (PSF) across the image 
by fitting each source as a 2-D gaussian. In Figure [6] we 
show the histogram of FWHM values. The high-end tail 
of the shown distribution is affected by truly extended 
objects, and a few odd noisy objects produce a roughly 
constant tail of objects with a wide range of FWHM. We 
use this tail to determine a non-zero offset, above which 
we fit the width of the remaining distribution as a = 0.13 
pixels. 

We assumed that all these sources are point-like stars 
and randomly placed them back onto an empty image 
with uniform noise consistent with the noise in our real 
image. We input stars as gaussians with the distribution 
of FWHM values derived above — a mean of 2.55 pixels 
and a standard deviation of 0.13 pixels. To test recovery 
of truly extended objects, we also inserted a relatively 
bright elliptical galaxy sliced out of our data frame ten 
times into each image. To simulate the additional selec- 
tion effects required to produce our catalog of sources we 
selected the 1024 objects with lowest photometric uncer- 
tainty (the same number of objects as in our real cat- 
alog) and retrieved their classification parameter value. 
The average number of objects assigned each classifica- 
tion parameter is shown in Figure [71 

We generated ten simulated fields and ran Source Ex- 
tractor with the same detection parameters used in the 
original image. In all cases pasted elliptical objects are 
recovered and properly classified as extended. However, 
we also found some other sources which are classified as 
extended. When we visually inspected these objects in 
the simulated fields, we confirmed that they were indeed 
superpositions of two points sources. The average num- 
ber of such objects in all simulations, is 4±2 or 0.4%. 
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Fig. 7. — The results of running Source Extractor on 10 test 
frames designed to reproduce the stellar density of Control 04. 
10 input galaxies per frame are recovered, while a small number 
of other objects are classified as extended (Class < 0.5). These 
objects are from overlapping point sources, and constitute roughly 
0.4% of the total objects. 

3.2. Extragalactic Candidates 
3.2.1. Quasars and AGN 

Active galactic nuclei could be classified as point-like or 
extended dep e nding on the brightness of the host galaxy. 
iHewett et al.1 ()2006l ) have produced color predictions for 
the WFCAM photometric system and conversions be- 
tween these and the 2MASS colors. From these colors it 
appears as if quasars with z < 6 might show up as near 
infrared excess sources. 

A careful calculation of the detecti on of quasars in the 
near i nfrared has been performed by iMaddox fc Hewettl 
(2006). We use their result for a if -band survey to 18.5, 
and assuming their most restrictive model for the com- 
bined quasar and galactic light (equivalent to being a 
point-like source), the estimate from their Table B5 is 
roughly 35 quasars per square degree up to 18th magni- 
tude, or roughly 2 in our 13 arc minute by 13 arc minute 
control field. 

3.2.2. Galaxies 

Surveys of galaxies in the near infrared are typically 
done in the K and J bands. The reason is that J — K 
has long been used as a simple galaxy-star discrimina- 
tor, using a criteri on where any object w ith J — K » 1 
is a galaxy (e.g. ISaracco et alj (|1999t )1. This distinc- 
tion corresponds to a diagonal line in our color-color di- 
agrams (see Figures &HJ:), and agrees marginally well 
with our structural classification. The 2MASS survey 
contains J, H and K data, but galaxies in that sur- 
vey are limited to objects much brighter than in our 
field. Few other surveys use H band , although this 
is be gi nning to chang e (c.f. UKIDSS (|Lawrence et al.1 
120061 ). iMetcalfe etall (|2006l )). Therefore, we consider 
number density estimates, a clustering analysis, and the- 
oretical k-corrections to understand the properties of our 
extended objects. 

Our density of extended objects agrees moderately 
well with de e p extr agalactic number density estimates. 
lElston et ail (|2006l ) provide data from the FLAMIN- 
GOS Extragalactic Survey (FLAMEX) and compare 
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number densities with previous surveys. Integrating 
over their differential number counts down to magni- 
tude 17.3 (our completeness of 17.4 with Ak = 0.07) 
gives roughly 3000 galaxies per square degree, or roughly 
140 galaxies detected in our field, compared to the 104 
extended objects we identify. A second estimate is 
avail able from the GALAXYCOUN TS 2 tool available 
from lEllis fc Bland-Hawthornl (|2007f ) which produces es- 
timated number counts and variances in several bands 
including K. This tool predicts 150 ± 27 galaxies for our 
survey's characteristics. This result is almost 2a dis- 
crepant with our numbers, but is quite sensitive to our 
completeness estimate and small variations in magnitude 
definitions. If we were only complete to a depth of 17.1 
then the estimate from lEUis fc Bland-Hawthorn! (|2007l ) 
of 121 ± 24 galaxies is la consistent with our number. 

In order to test that our control field was not point- 
ing by coincidence at a cluster of galaxies, we compared 
the surface density of extended objects (C < 0.5 and 
K < 17 A) wi th that of galaxies in one of the FLAMEX 
survey fields (|Elston et a l. 2006) . We selected from their 
Cetus catalog one region named R3D3 with central co- 
ordinates (a(2000),<5(2000))=(33. 91840,-4.67243) which 
is known to have a prominent galaxy cluster; we se- 
lected galaxies with K s < 17 .4 and ellipticities (1 — b/a) 
larger than 0.2, in order to avoid including any signif- 
icant numbers of stellar profile objects, which resulted 
in a total of 305 galaxies. To compare the surface den- 
sities, we calculated 6th nearest neighbor distances, de, 
for objects in both lists, which c an be translated to local 
surface densities as 5/ (ird%) (c.i. ICasertano &: HutlH985t 
iFerreira fc Ladall2007f ) . In the Cetus-R3D3 and the Con- 
trol 04 fields the mean densities of objects are 1.06 ± 1.30 
and 0.53±0.25 objects per square arcmin, while the max- 
imum densities are 11.14 and 1.33 objects per square ar- 
cmin, respectively. In Figure [5] we show maps of 6th 
neighbor surface densities for these fields in steps of 0.5 
objects per square arcmin. While the nearest neighbor 
map clearly identifies the Cetus R3D3 cluster, the map 
for Control 04 shows no evidence of a cluster. 

Finally, we seek to understand why our (non-cluster) 
galaxies form a relatively tight clump in the near-infrared 
color-color diagram. Down to if =17.4 mag, the depth of 
our observations, a significant number of galaxies would 
be expected to be at moderately high redshift. Two ef- 
fects can contribute to change the colors of distant galax- 
ies. The first is the fc-correction, which means that at 
higher redshift our near-infrared bands are observing in- 
trinsically bluer portions of the galactic spectrum. The 
second is intrinsic galactic evolution — different star for- 
mation rates, for instance. 

Due to the paucity of TJ-band deep extragalactic imag- 
ing, we do not have a way to compare colors of genuine 
red-shifted galaxies with those of the objects in our NIRX 
clump. Thus, we rely on theoretical fc-corrections to un- 
derstand the colors of our extended obje cts. We consider 
two di fferent k-correction models, that of lMannucci et all 
(|200l[) based on near-inf rared template spectra, and the 
HYPERZ templates by iBolzonella et alj (|200dD . which 
are visual spectra extended in wavelength using a popu- 
la tion synthesis code. Both these models are compared 
in lHewett et al.1 (|2006l ). who also carefully consider their 

2 http:/ /www. aao.gov.au/astro/GalaxyCount/ 



transformations between different filter systems, and we 
use their transformations into the 2MASS photomet- 
ric system 3 . Also, we apply a constant extinction of 
Ak =0.07 to both models to account for the extinction 
present in our control field. 

The comparison between our data and the fc-correction 
models is shown in Figure [51 We found that there are 
fairly significant differences in the shapes of the two mod- 
els, but most of our extended objects can be explained as 
galaxies at moderate redshift. To quantify this, we bin 
our extended objects by if-band magnitude into three 
magnitude bins from 14 to 17. For each bin, we mea- 
sure t he median and la dis persion of H — K and J — H. 
From ISongaila et ail (|1994h we obtain the expected me- 
dian redshift for each if-band magnitude bin. This red- 
shift is converted into an H — K and J — H estimate 
using both the Kinney-Mannuci models and the HY- 
PERZ models and the results shown in Table [2] The 
HYPERZ templates do not match our observed colors 
(missing most severely in the suspect J band), but the 
Kinney-Mannucci templates appear to explain well the 
colors of our NIRX sources as being due to a distribution 
of normal galaxies of various types at redshifts between 
0.1 and 0.5. 

We are not able to distinguish this population of nor- 
mal galaxies from more exotic types of AGN-dominated 
or starburst galaxies based on NIR colors alone. How- 
ever, since our number density estimate for normal galax- 
ies is larger than or consistent with our observed number 
density, and we understand the color-color clumping be- 
havior as a consequence of the fc-correction, we conclude 
that we are not seeing a significant population of exotic 
galaxies. 

4. IMPLICATIONS FOR CLUSTER STUDIES 
4.1. Cluster Studies in the Near-infrared 

Our analysis has a direct implication in the study of 
star forming regio ns. We consider data from the Rosette 
Molecular Cloud ( Roma n- Zuniga et al., 2007) which is of 
insufficient quality (due to instrumental problems) to 
make a good distinction between extended and point- 
like objects. However, it illustrates how contamination 
by galaxies could be important in cluster studies. 

J — H vs. H — K color-color diagrams are a common 
tool used to determine which objects in a certain region 
have circumstellar emission. In the classical picture, the 
data are plotted in the color-color diagram and every ob- 
ject falling to the right of the reddening band of the zero 
age main sequence (e.g. J — H < 1.692(H— K)+B, with 
B accounting for a typical color error buffer) is consid- 
ered a candidate young star. However, our study shows 
that a significant fraction of these objects, especially faint 
ones, may be galaxies, rather than YSOs. 

In many cases the quality of the observations can add 
to the problem. In order to make a good separation of 
elongated or fuzzy objects which might be extragalactic 
contaminants in a field populated with young stars us- 
ing automated software, seeing and resolution need to 
be good. Still, in some cases, even with good quality 

3 Our J-band magnitudes are taken with a filter which is sig- 
nificantly different from the 2MASS one. Therefore, although we 
calibrate to the 2MASS system, this transformation is inexact for 
non-stellar objects 
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TABLE 2 

Comparison of galaxy colors with models. 



Data properties HYPERZ model Kinney-Mannuci model 



Mags 


Median(z) 






Median 


17 


Elliptical 


Sbc 


Scd 


Ellip. 


Sa 


Sb 


Sc 


14-15 


0.18 


H - 


K 


0.53 


0.10 


0.627 


0.588 


0.570 


0.543 


0.570 


0.583 


0.562 






J - 


H 


0.82 


0.10 


0.701 


0.672 


0.647 


0.881 


0.865 


0.897 


0.87 


15-16 


0.28 


H - 


K 


0.70 


0.14 


0.739 


0.688 


0.665 


0.706 


0.720 


0.704 


0.701 






,1 - 


H 


0.90 


0.09 


0.766 


0.732 


0.668 


0.938 


0.905 


0.941 


0.896 


16-17 


0.35 


H - 


K 


0.73 


0.16 


0.729 


0.688 


0.665 


0.756 


0.765 


0.747 


0.748 






J - 


H 


0.95 


0.14 


0.835 


0.801 


0.692 


0.962 


0.931 


0.938 


0.891 



Note. — On the right are observed median col ors within set magnitude bins. By using the measured 
magnitude-redshift relation of Songaila ct al. ( 1994) we predict colors using the two different sets of models 
introduced in the text for various galaxy types. 
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Fig. 10. — Near-infrared color -color diagram for a contro l 
field in the Rosette nebula from Roman-Zuniga ct al. (2007), 
with the same lines and regions as Figure [T] Unable to use 
structural information to e liminate the clump of NIRX sources, 
I Roman -Zuniga ct al. (2007) were forced to trim at K = 15.75 mag- 
nitudes (points fainter than this limit are shown as grey stars), 
roughly 2 magnitudes brighter than their detection limit. This 
tightens the dispersion on stellar colors and eliminates most of the 
NIRX sources. 



observations, some extragalactic objects might not be 
detectable if their elongations are small, and thus some 
sources may pass as genuine young stars. The best so- 
lution, in most cases, is to have a good quality control 
field near the region of interest and perform in it a simi- 
lar analysis to the one we performed in the Perseus field, 
which can give an estimate of the level of contamination 
in the field. Other factors like extinction might also have 
to be taken into account especially if the star forming 
region is very young. 

In their study o f the Rosette Molecular Cloud, 
iRoman-Zuniga et al.l ((2007) decided that for sources 
fainter than K = 15.75 mag they could not distinguish 
genuine NIRX sources from galactic or extragalactic con- 
taminants, or from distorted objects produced by instru- 
mental problems. Shown in Figure [10] is one of their 
control fields displayed as in Figure [1] showing both faint 
(grey) and bright (black) sources. A similar set of NIRX 
sources is seen, mostly in faint objects. However, us- 
ing a slightly larger area for possible NIRX sources than 
we use in this paper, they found that even after cutting 
out faint objects their two control fields (carefully se- 
lected to be located outside of the Rosette Cloud) still 
contained some NIRX sources They used this number of 
NIRX sources to define a background level, which helped 
them to determine which regions of the Rosette have 
young cluster formation (regions with surface densities 
of NIRX sources above the defined background). If they 
had not corrected by this background population, then 
larger area of their cloud fields would present a high den- 
sity of infrared excess sources and might have passed as 
star forming regions. 



4.2. Identifying YSOs with Spitzer data and 
High-quality NIR Data 

Galaxies masquerading as YSOs have been a signif- 
icant headache for the analysis of Spitzer observations 
of star-forming regions, and contamination by galaxies 
may set the limit for identifying the lowest luminos - 
ity YSOs from these data alone (|Harvev et al.l [2007) . 
Near-infrared imaging is but one source of additional 
data which may be brought to bear on this problem. 
We combine our observations with data from the Cores 
to Disks (c2d) Legacy program survey of Perseus pre- 
sented in |j0rgensen et all ([2006) to check our classifica- 
tion scheme and to demonstrate the value of adding high- 
quality near-infrared imaging to these rich datasets. As 
we will show, additional constraints come mainly from 
structural, rather than color information. Unfortunately, 
the refined distinction between point-like and extended 
objects provided by NIR data is still not a definitive way 
to separate YSOs from galaxies, as YSOs may exhibit 
outflows or nebulosity, and distant galaxies may remain 
unresolved. 

No Spitzer data e xist for our control fields, but 
lJ0rgensen et al.l (|2006( ) present c2d program data for all 
six on-cloud frames. Following their procedure for iden- 
tifying YSOs, we cross-correlated our NIR data with the 
c2d catalog (DR3) and pulled out magnitudes for ob- 
jects which were point-like with photometric quality A 
or B in all four IRAC bands. If a high-quality detec- 
tion was made at 24 /zm(MIPS), that information was 
also included. B ased on a compar i son w ith a SWIRE ex- 
tragalactic field, lJ0rgensen et al.1 (|2006l ) give two differ- 
ent color-magnitude cuts to identify YSOs. Points lying 
above and to the right of these cuts in either the [8.0] 
vs [4.5]-[8.0] or the [24] vs [8.0]-[24] color-magnitude di- 
agram are considered candidate YSOs. These cuts are 
shown in Figure 1111 along with our classification of ob- 
jects as extended (open symbols) or point-like (filled sym- 
bols). The well-known Class protostars in L1448 and 
B5 are excluded from this catalog as they have outflows 
making them appear extended in the IRAC images. 

IRAC images are limited to a pixel scale of 1.2" and 
roughly 2" resolution. With our 0.45" pixel scale and 
0.6" seeing, we are able to better identify marginally 
extended objects than Spitzer, as every object in Fig- 
ure[Tl]is unresolved by Spitzer. The structural and color- 
magnitude divisions of objects into stellar and extended 
(galactic) objects matches exceedingly well for most of 
the 402 objects, but there are ten interesting objects. 
The EE gensen et al.l (|2006h criteria identify 5 candidate 
YSOs. Of these, three are point-like and bright; these 
are most likely genuine YSOs (purple diamonds in Fig- 
ure [TTj) . Two candidate YSOs are significantly extended 
and thus more likely to be galaxies (green crossed-boxes 
in Figure ITTj) . Five objects unresolved even in the NIR 
have galaxy-like colors and magnitudes (green asterisks 
in Figure ITTj) . Of these, one is relatively extended, and 
the other four are all much fainter than the other point- 
like sources. We cannot definitively establish the nature 
of these objects — they could either be unresolved galax- 
ie s or low-luminos i ty YS Os . 

|j0rgensen et al.l (|2006h combined 2MASS with Spitzer 
data in both their Perseus data and an extragalactic 
survey and based on this suggest that a H — K vs. 
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Fig. 11. — Based on NIR data for six on-cloud frames we show (a) structural information for objects with high-quality detections 
in all four bands in IR AC which appear point-like; we resolve many as extended. Of the five objects identified as YSO candidates by 
|j0rgensen et all (2006) based on either color-magnitude diagram (b) or (c), three are point-like and likely YSOs, and two are extended and 
probably galaxies. A handful of potentially point-like objects are either unresolved galaxies or low-luminosity YSOs previously unidentified. 
The addition of H and K photometry (d) is ins ufficient to distinguish between YSOs and galaxies. The shaded region is the area where 
"only YSOs" are found in j0rgcnsen ct al. (2006) using 2MASS data, but this is because 2MASS only picks up bright local galaxies. With 
deeper near-infrared data, the "YSO-only" area is littered with galaxies. 



K — [4.5] color-color diagram should provide a clean sep- 
aration between YSOs and galaxies. Unfortunately, this 
appears to be due to the magnitude limit of 2MASS, 
which only picks up local galaxies with star-like colors. 
With deeper NIR data, the region of the color-color dia- 
gram where Spitzer and 2MASS see only YSOs is actually 
heavily contaminated with galaxies (Figure [TTtO. Struc- 
tural, rather than photometric, information appears to 
be be required for making the YSO/galaxy distinction 
for sources fainter than K — 15. 

5. IMPLICATIONS FOR EXTINCTION MAPPING:GNICER 

In studies of young clusters near-infrared excess sources 
constitute both signal (YSOs) and noise (everything 
else), but in extinction mapping using background stars 



all NIRX are potential contaminants. Brown dwarfs are, 
by and large, foreground sources, with colors unrelated 
to the cloud material. YSOs may be associated with 
ongoing star formation and appear in front of, within, 
or behind the cloud tracing none, some, or all of the 
cloud material. Other near-infrared excess sources are 
reddened by the cloud, but since their intrinsic color is 
redder than normal stars they will invariably produce an 
overestimate of the column density. 

However, as we have shown above, in regions relatively 
far from the galactic plane the majority of near infrared 
excess sources are galaxies, and therefore they are behind 
the cloud and potentially useful tracers of column den- 
sity As both the noise and resolution of extinction maps 
depend on the number of background sources, finding 
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TABLE 3 

Comparison of Different Extinction Algorithms 



No. 



Algorithm 



Shown in Figures... 



Discarding galaxies (NICER) 
Treating all sources as stars 
Using different mean colors and dis- 
persions for galaxies and stars 
Same as (iii), but using magnitude- 
dependent properties for galaxies 
(GNICER) 



14,15,16 (x-axis) 17 
14 (y-axis) 17 
15a (y-axis) 

15b,16 (y-axis) 17 



a way to use these galaxies is desirable. Indeed, in one 
way galaxies are an id eal probe, since they a re never fore- 
ground to our cloud. lYasuda et al.l (|2007f ) has recently 
used ga laxy number coun ts from the Sloan Digital Sky 
Survey ijYork et all l2"f)0l as a probe of extinction. We 
seek to make use of galaxy NIR colors, which we expect 
to be a more sensitive measure. 

We consider and compare four approaches to making 
an extinction map of L 1451, a starless cloud in the south- 
west of Perseus: i) eliminating galaxies, ii) the naive 
treatment of blindly applying near-infrared extinction 
mapping to all sources, and iii-iv) two different methods 
for using galaxies. These different algorithms are sum- 
marized in Tabled The color-color diagram for L1451 is 
shown in Figure [12] and in all cases we use Control 04 to 
determine the intrin sic properties. A JHK col or image 
of L1451 is shown in lFoster fc Goodman! (|2006f ). 

As shown in Figure [T2l it is not always easy to identify 
galaxies in a reddened field based simply on color. If wc 
had J band detections for every object, we could make 
a cut in color-color space parallel to the reddening vec- 
tor and thus identify NIRX sources. However, reddening 
hopelessly mingles stars and galaxies if one only has their 
H—K colors, as is often the case behind high column den- 
sity regions. It is also possible to identify NIRX sources 
by using a preliminary extinction map to construct an 
intrinsic color-color d iagram of all sources. T his was the 
approach adopted bv lLombardi et aTl (|2006f ). They con- 
structed an initial extinction map for all their sources, 
converted observed colors into intrinsic, and eliminated 
near infrared excess sources. In our case, we are able 
to make a distinction based on structural information, 
with the "sigma clipping" algorithm described below re- 
moving the few objects which are extended with normal 
colors, or point-like with near infrared excess. 

Our baseline for comparison is an extinction map 
made with the standard application of the Near-Infrared 
Color Excess Method Revisited (NICER) after elimi- 
nating galaxies (Method i). Following the notation of 
lLombardi fc Alvesl ()2001l ). each object provides a pencil 
beam measure of the extinction from 



A v = bj-„[(J-H)° bs -(J-H) con } 
+ b H _ K [(H - K) obs -(H- K) con 



(2) 



Assuming a normal reddening law, bj-u and bu-K are 
constrained such that 



bj-H.f,j-H + ^h-kIh-k — 1 



(3) 



where Jj-h is the ratio of color excess to visual ex- 
tinction, fj-H = E( J — H)/Av = 1/9.35 and likewise 
Sh-k = 1/15.87. 
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Fig. 12. — L1451, a cloud with no known embedded young stars 
in the southwest of Perseus. Red open square are extended and 
blue closed symbols are point-like. Point-like and extended objects 
without a high quality J detection are plotted at J-H of and 0.1 
respectively to show their distribution in H-K. A fit to the red- 
dening vector (shown) is consistent with that offRickc & L cbofskvl 
1119851) . 

As the colors J—H and H—K are not independent (i.e. 
there is a stellar locus in the color-color diagram) and 
the photometric errors on the colors are also correlated 
(though we assume that errors on individual filters are 
uncorrelated) we calculate the covariance matrix C — 
E + G, with 



E 



' H 



H 



' H 



H 



K 



(4) 



the covariance of the photometric errors and 



r _( Var( J - H) 

{ Cov( J - H,H — K) 



Cov( J ~H,H-K) 
Var(H - K) 



( 5 ) 

the intrinsic relationship in object colors measured from a 
control field. We wish to minimize the variance of Eqn.[3 
subj ect to the constraint of E qn. [3] This is equivalent to 
(see [Lombar di fc Alvesl (|200lD for full details) solving 



C" 1 f 
f -C- 1 i 



(6) 



with b = (bj- H ,b H ^ K ) and f = (fj-H,fH-i<)- The 
error associated with each point like determination is 



Var(Av) = b • C • b. 



(7) 



A map is created by smoothing these point estimates 
with a gaussian beam, weighting the contribution to a 
cell's extinction determination by the error on the ob- 
ject. Thus the error on a cell's extinction value is roughly 
inversely proportional to the square root of the number 
of point estimates used in that cell's determination. In 
the densest portions of the cloud, a significant number of 
background objects are not observable, and so contami- 
nation by foreground stars becomes a problem (typically 
this becomes a problem around Ay = 10 mag (|Lombardil 
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Fig. 13. — A point-to-point comparison of an extinction map for 
L1451 derived naively assuming that every object has the mean 
color of the control field (Method ii), versus just using stars (i.e. 
point-like) objects (Method i). The solid line shows a one-to-one 
relation, and the dashed line shows a fit to the relation. 

2005)). To reduce this bias we perform "sigma clipping" . 
To do this, we determine a pixel's extinction, clip points 
more than three sigma deviant from the estimated Ay, 
and recalculate Ay until this process converges. Where 
fewer than 5 objects contribute to a pixel's extinction 
determination we assign a blank value in the map and 
exclude that pixel from further analysis. 

By rejecting extended objects in both our data frame 
and our control field property determination, we create 
our standard map for comparison (Method i) . Figure [T3] 
shows the comparison between this map and the naive 
approach in which NIRX sources are not accounted for 
(Method ii), and so all sources are considered to be stars 
and de-reddened to the stellar locus. This naive approach 
is biased towards determining higher extinctions, as in- 
trinsically red galaxies give large estimates of the column 
density. The problem is mitigated at high extinction be- 
cause the galaxy population is intrinsically fainter than 
the stellar one and thus becomes a decreasing percent- 
age of the sources seen through a thick cloud. Both these 
maps have a pixel size of 45" and no holes (i.e. pixels 
with less than 5 objects contributing). 

Since we can separate stars and galaxies without re- 
course to their colors, galaxies provide crucial extra in- 
formation in making an extinction map. We therefore in- 
troduce the Galaxies Near Infrared Color Excess method 
Revisited: GNICER. The basic approach (Method iii) is 
to divide objects into point-like and extended classes as 
before in both data and control fields and determine pa- 
rameters for these populations separately. That is, solve 
for (J — H) con and (H — K) con for use in Eqn.[2]for both 
point-like and extended objects, and also determine the 
intrinsic covariance matrix G in Eqn. [5] for both popula- 
tions separately. Every source now gives a point-beam es- 
timate with the uncertainty on that estimate (Eqn.[7| re- 
flecting both the photometric error and how tightly such 
objects are intrinsically clumped in color-color space. Be- 
cause our galaxies are approximately as tightly clumped 
as our stars, each galaxy provides roughly as much infor- 
mation about extinction as a star (subject to photometric 



uncertainty) . 

The most stringent requirement is that an extinction 
map created from only galaxies should provide roughly 
the same result as one created from just stars. The re- 
sult of comparing GNICER (Method iii) for just galax- 
ies with the straightforward application of NICER (for 
point-like sources only, Method i) is shown in Figure[T4"h. 
Note that especially at larger Ay the galactic estimate 
becomes quite scattered as there are fewer galaxies. As 
before, this analysis excludes pixels in the extinction map 
with fewer than five objects contributing to the estimate 
within that pixel (this is a number of pixels in the galaxy- 
map). The linear best-fit to the data (including errors in 
both X and Y) reveals a slope which is statistically sig- 
nificantly less than unity. Since the galaxies are detected 
only down to a brighter limit, this bias could arise from 
extinction structure within the pixel — the fainter galax- 
ies trace low density material compared to the brighter 
stars in our magnitude limited image. Another source 
of this bias is the fact that brighter galaxies are pref- 
erentially bluer as fainter galaxies are typically further 
away, have a larger fc-correction, and are redder. Behind 
significant extinction, we see only these bright, intrinsi- 
cally bluer objects. Therefore, we estimate that they are 
sitting behind less extinction than they really arc. 

We test this second source of bias in our fourth exper- 
iment (Method iv) by splitting the control field galaxies 
into magnitude bins and determining the mean colors 
and color covariance matrix for galaxies up to a certain 
magnitude, Ki 00 k U p- We make a rough first extinction 
map using just stars. Then, in each pixel with extinction 
of Ay, galaxies are used to refine the extinction deter- 
mination. For a galaxy of magnitude K i, s we determine 
which set of galaxy properties from the control field to 
use by using 

Ki 00 kup — K obs — 0.112Ay. (8) 

The properties of our galaxies binned this way is shown 
in Table 3J Figure [14b show the result of this more com- 
plex algorithm, again comparing an extinction map made 
just with galaxies (y-axis, Method iv) against one made 
just with stars (x-axis, Method i). Most of the bias 
(non-unity slope) is removed by considering the color- 
magnitude relation of the galaxies. 

When we combine the galaxy and stellar estimates in 
GNICER (Method iv) we make a map which is unbi- 
ased with respect to the map which eliminates galaxies 
(Method i, Figure [T5"]) . The big win is in the noise of the 
map, which is significantly reduced by adding galaxies as 
crucial extra estimators of the extinction. A comparison 
of the error distribution between naively assuming every- 
thing is a star (Method ii), rejecting galaxies (Method 
i), and this second algorithm for including galaxies in 
GNICER (Method iv) is shown in Figure Q3J The naive 
treatment has relatively low estimated errors since it too 
uses all the background sources, but this is misleading 
because it is biased - estimating much higher extinction 
for regions with low intrinsic extinction. 

Future work which attempts to effectively use galax- 
ies as additional sources for extinction mapping will ide- 
ally be able to: a) confidently categorize sources based 
on structural information (this requires adequate seeing, 
pixel scale, and image quality) , and b) have control fields 
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Fig. 14. — A point-to-point comparison of an extinction map for L1451 as derived from just galaxies versus just stars (Method i). 
Two different algorithms are compared. The simple treatment in (a) treats all galaxies as having the intrinsic properties of the control 
field galaxies (Method iii). The complex treatment in (b) assigns an intrinsic color to a galaxy based on its intrinsic magnitude using an 
extinction estimate from the stars (Method iv). This complex treatment helps remove bias at high Ay. The solid line shows a one-to-one 
relation, and the dashed line shows the fit. 



TABLE 4 

GN1CER PARAMETERS 



Depth 


(J-H) 


(H-K) 




Gl2 


G22 


Stars 












All 


0.584 


0.177 


0.027 


0.010 


0.017 


Galaxies 












15.0 


0.817 


0.530 


0.010 


0.010 


0.010 


15.5 


0.856 


0.635 


0.006 


0.006 


0.016 


16.0 


0.878 


0.671 


0.011 


0.012 


0.026 


16.5 


0.914 


0.695 


0.024 


0.021 


0.028 


17.0 


0.935 


0.696 


0.021 


0.016 


0.027 


17.5 


0.933 


0.671 


0.019 


0.011 


0.030 



Note. — Measured mean colors and color-covariance 
matrices of stars and galaxies binned by magnitude 
from Control 04. These parameters are used in 
GNICER to incorporate galaxies into an extinction 
map. Gn = Var(J - H), G12 = Cov(J - H,H - 
K),G 2 2 = Var(H - K) 

with adequate depth to determine the color-magnitude 
dependence of galaxies as seen through whatever portion 
of the galaxy is being studied. Galaxy clusters, and large- 
scale structure in the cosmological sense, pose a problem 
for GNICER. A galaxy cluster will be a dense group of 
objects at a small redshift range, perhaps significantly 
different from the average redshift for all galaxies. How- 
ever, it is not clear that this problem is more severe than 
the problem which stellar clusters present NICER. Spa- 
tially coincident clumps of abnormally colored objects 
must be identified and understood in either case. 

6. CONCLUSION 

We have identified a number of near infrared excess 
sources (NIRX) in control fields in Perseus as normal 
galaxies at moderate redshifts. With high-quality im- 
ages (i.e. good seeing conditions and fine pixel scale) we 
are able to classify objects as stars or galaxies based pri- 
marily on structural information. In such fields, galaxies 
which appear as NIRX sources constitute a significant 
background level in surveys for YSOs, and a significant 
potential contaminant in extinction maps using NICER. 




10 

A v using just stars 

Fig. 15. — A point-to-point comparison of an extinction map 
for L1451 derived by separating objects into two classes (extended 
and point-like) and using the magnitude-binned properties (mean 
color and dispersion) of these two populations as determined from 
the control field (Me tho d iv), versus using just stars (Method i). 
In contrast to Figure fill the y-axis shows points from a map made 
from both stars and galaxies, rather than just galaxies. The solid 
line shows a one-to-one relation, and dashed line shows the fit. 



High-quality NIR imaging of cloud complexes observed 
in the mid-IR (e.g. with Spitzer) should provide the nec- 
essary discriminant to separate YSOs from galaxies in 
these data sets, extending the YSO luminosity function 
at the faint end where galaxy contamination currently 
dominates. This separation must be made from struc- 
tural information, rather than relying on only colors. 

We have extended NICER to include galaxies in an un- 
biased way despite galaxies' significant color-magnitude 
relation (GNICER). Including galaxies in the analysis of- 
fers a significant increase in the number of background 
sources used to probe the structure of the cloud and a 
corresponding decrease in the uncertainty of this esti- 
mate. Using GNICER it should be possible to map low 
density clouds relatively far from the galactic plane with 
reasonable resolution. 
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Fig. 16. — A comparison of the map error amongst the three 
different way of making an extinction map (naively treating every- 
thing as a star, throwing out galaxies, and using galaxies). Errors 
on the naive determination (Method ii) are low, but this map is 
biased with a non-zero intercept and non-unity slope in compari- 
son with a map where galaxies are discarded (Figure 1131 Method 
i). Including galaxies (GNICER, Method iv) gives many additional 
sources from which to estimate the extinction. Thi s reduces the 
error without introducing significant bias (Figuro ll5H . 
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